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Abstract 



A doublet of three-dimensional Dirac fermions can effectively describe the low en- 
ergy spectrum of a fermionic cubic lattice. We employ this fermion doubling to 

Cd ' encode a non-Abelian SU{2) charge in the fundamental representation. We explic- 

^ . itly demonstrate that suitable distortion of the tunnelling couplings can introduce a 

i-rt ! scalar and a Yang- Mills field in the effective low energy description, both coupled to 

Ch ' the Dirac fermions. The simplicity of the model suggests its physical implementation 

^ ■ with ultra-cold atoms or molecules. 
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. . 1 Introduction 
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oo 

O ■ Recently, much interest lias been focused on fermionic lattice models that can 

simulate one or two-dimensional exotic phenomena such as fractionalization 

^ ', of charge or anyonic statistics. Polyacetylene is a one-dimensional example, 

H I which gave a simple and rigorous theoretical model for the charge fractionaliza- 

tion [Ij that enjoyed experimental verification [2J. In two-dimensions graphene, 
a mono-atomic layer of graphite, has a low energy description given by a two- 
dimensional Dirac equation [3J. It has been shown that under certain coupling 
distortions of the underlining lattice it can support axial gauge theories that 
lead to the fractionalization of charge [H[5|S|7] . In three or four dimensions the 
properties of these lattices and their coupling to external non-Abelian gauge 
fields are well understood, for example, in the context of staggered fermion 
lattice gauge theories [T0llllfl2j . How non-Abelian gauge fields can effectively 
arise in the continuum limit without employing them already in the lattice 
picture remained up to now an open question 
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In this letter, we systematically study the case of three spatial dimensions. 
We employ a cubic lattice model where fermions are tunnelling from site to 
site with specific transition couplings. It has been shown [13j that the low 
energy spectrum of this model is effectively described by three-dimensional 
Dirac spinors. Here, we demonstrate that suitable distortions of the transition 
couphngs and of the local potentials give rise to additional triplet scalar and 
Yang-Mills fields, both of them coupled to the Dirac field. Importantly, the 
presented distortions do not originate from non-Abelian fields in the lattice 
picture. The interest in such structures is diverse. The simplicity of these 
lattice models could inspire physical systems such as optical lattices or multi- 
layer graphene that can simulate Dirac fermions coupled to dynamical Yang- 
Mills gauge fields. Such an endeavor could establish a bridge between the 
lattice gauge theory and the experimental condensed matter communities. 
Moreover, it is of interest to see if the standard model can emerge as an 
effective theory of a simple underlined lattice structure [9j . Such a model could 
generate Abelian and non-Abelian gauge potentials as well as Higgs scalar 
fields that emerge without a symmetry breaking due to a simple quantum 
ordered in its ground state. 



2 The three dimensional case 



2.1 The fermionic lattice model and the Dirac equation 



For concreteness we consider the three-dimensional Hamiltonian 



D^ = (l®7-p-T'^®7-A, + T'^®l<l>,)^, (1) 

where summations are assumed over repeated indices. The gauge field is taken 
to be in the fundamental representation of SU{2), i.e. T"" = ^ [14j. The four 

spinor, ^^ , is indexed by z = 1,...,4, while n = 1,2 runs trough the color 
components. The index a = 1,2,3 runs through the gauge field components 
and {7^,7;/} = 26^^^, where fi, v run through the space components x, y and 
z. By choice of gauge the three-dimensional gauge vector A^ and the triplet 
scalar field $« are taken to be real. 

Our starting point is a cubic lattice where fermions tunnel from one site to 
the neighboring one according to the simple quadratic Hamiltonian 

i^ = -^V^|XijV^j + H.c.. (2) 

<iJ> 
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Fig. 1. A two-dimensional section of the three-dimensional cubic lattice. Fermions, 
V'i, tunnel from site i to site j along the links of the lattice with couplings 
"^ij ~ Xij ^ "^Xij + ^Xij- Here Xij and ^Xij couple nearest neighboring sites, i and j, 
that lie only along the same link of the cubic lattice, while Axij couples diagonal 
sites as well. 

The tunnelling couplings are given by 



^iJ 



Xij + ^Xij + Axij, 



(3) 



where Xij correspond to a uniform pattern, while (5xij and Axij correspond to 
coupling distortions that could have a slow spatial variation. While Xij and (5xij 
connect nearest neighbor sites that lie only along the links of the cubic lattice, 
the couplings Axij can also generate site energy shifts or couple sites along 
diagonals, as seen in Fig. [H As we shall see the couplings Xij are chosen so that 
the low energy description of H is given by a Dirac equation. There are three 
independent components of (5xij, along the different directions of the lattice, 
that can give rise to a triplet scalar field, while nine independent components 
of Axij can be identified that can give rise to the SU{2) non-Abelian gauge 
potential. 

Let us first consider the effect of the uniform couplings Xij that give rise to 
the Kogut-Susskind (or staggered) fermions [13]. Their particular values that 
actually minimize the energy of Hamiltonian ([2]) are given by 



Xi,i+S 



:) Xi,i+y 



;-l)*^ xi,i+z 



'_^Yx+ly 



This configuration corresponds to a vr-flux going through each square plaque- 
tte of the lattice. As the couplings, Xij? have alternating signs along the x and 
y directions, while they are homogeneous along the z direction, the result- 
ing Brillouin zone, B, is given by kx,ky G [— 7r/2,7r/2] and k^ G [— 7r,7r]. In 
momentum space the Hamiltonian kernel, H, is given by 

iJ = 1 (g) (jg sin fcj. + (73 CTi sin ky + cxi ® o"i sin k^ (4) 

in the basis (k, k + Qx, k + Qy, k + Qx + Qy), where Qx = (tt, 0, 0) and Qy = 

(0, vr, 0). The corresponding eigenvalues of the energy are given by 



E, 



±\/sin kx + sin kj, + sin kz- 



Within the Brillouin zone, B, there are two isolated Fermi points given by 
K_|_ = (0,0,0) and K_ = (0,0, vr) that set the energy, E±, to zero. Focusing 
on the low energy regime of the Hamiltonian at half filling we expand (jlj) 
around K_|_ and K_. For k = K-t + p we have 

H± = 1 ^ a3p^ + (T3 ® aiPy ± (Ti ® (Tip^, (5) 

where |p| is taken to be small. We can now define 7^: = 1 (g) o"3, •jy = as ® ai 
and 72 = (Ti 0"! so that Hj^ = 7 ■ P- The corresponding spinors are given 

by ^±(p)^ = (V^±,p,V^±,p+q.,V'±,p+q,,^±,p+q.+qJ where V'i.p = V^(K± + p). 
Since ((73 ® l)i/_((73 ® 1) = iJ^ it is easy to show that for U = diag(l 
l,cr3 




= 1®7-P- (6) 

In this way one can obtain two identical copies of the free three-dimensional 

Dirac equation. The spinor, \E', of eqn. ([1]) is given hj ^\ = f/^E'+j and 

(2) 
\E'j- = t/^E'-j, where i = 1, ..., 4 enumerates the four components of the spinors. 



2.2 The scalar field 



In order to obtain the scalar field of Hamiltonian ([T]) we choose the coupling 
distortions to be of the form 

Sm+y = -|e-^^-'e^^+-('+^^), 5xi,i+z = -y , (7) 

where $ = $1 — i$2 with all $„ real. The i vector indicates the position of the 
corresponding site and the s^ vector indicates its two nearest neighbors along 
the direction a. The components of the scalar field can have a variation in 
space, which is slow compared to the lattice spacing. As 5xi,i+x = we have 
that SH^^^ = 0. For the distortion Sxi,i+y and for momenta close to the Fermi 
points we obtain the additional Hamiltonian term 

<^tt(ij) $10"! + $20-2 

The distortion in the z direction, 5xi,i+z, gives the Hamiltonian contribution 

Finally, the conjugation U{SH^y^ + SH^^^)W results to the triplet scalar field 
term, T" ® l$a, of Hamiltonian ([T|). 



2.3 The Yang-Mills field 



The gauge field can originate from tlie Axij distortion. Tlie non-zero couplings 
are given by 



A^ Ay 

AX 

^ -JK_i„iK+i 



Axi,i+x+y = -— e '"^-'e 



1 



A^ A^ 

Axi,i+x+z = -y|, Axi,i+y+z = -(-1)'^^' (8) 

where A^ = Ai — iA2 with all A'^ real and slowly varying with respect to the 
lattice spacing. These distortions of the tunnelling couplings provide the nine 
components of the non-Abelian gauge potential. Indeed, Axi i+y gives rise to 

From Axi i, which is a distortion of the local potential of site i, we obtain 

From Axi i+x4.^ we obtain the additional Hamiltonian term 

while from l^X\w±^i we obtain 

^^(x+z) _ _^^^ ® 1 ® as. 
Finally, from Axi,i+y+z we obtain 

After conjugating by \J we have that U{AH^y^ + A/f (°) + AiJ(^+y) + Aif (*+^) + 
Aif(y+^))f/t = —2"" (g) J . Afl. This shows that the Axij coupling distortions 
followed by the U conjugation give rise to the gauge potential term of Hamil- 
tonian ([1]). Thus, by appropriate distortions of the tunnelling couplings of the 
fermionic lattice model one can generate a non-Abelian gauge and a scalar 
field coupled to a Dirac field. Note, that, while these distortions may have a 



complex value the total additional plaquette flux they introduce is zero. Thus, 
the initial lattice Hamiltonian ([2]) is time symmetric. 



2.4 Yang-Mills monopoles and the index theorem 



In our construction the gauge and scalar fleld components can take arbi- 
trary conflgurations. As an application, we would like to simulate Yang-Mills 
monopoles in the low energy behavior of the presented fermionic lattice model 
by appropriately choosing the coupling distortions, ^Xij a.nd Axij. The fleld 
conflgurations, induced by a single monopole positioned at the origin with 
unit charge, are given by 

$, = faHr), A^: = e^^^r.Air). (9) 

Both $(r) and A{r) are taken to vanish at r = 0, while for large r the scalar 

fleld tends to its vacuum value, $ > m, and the gauge fleld behaves as 

A > —I/'", where the approach to their asymptotic behavior is exponen- 
tial [15]. Substituting these conflgurations into Hamiltonian ([T]) one obtains 



D-^ 



l®7-p-y ®(7xr)„v4+y ®lf„$]^, (10) 



which describes the Yang-Mills monopole coupled to a Fermi fleld. 

The appearance of zero fermionic modes is the most striking characteristic 
of such non-trivial topological conflgurations [16]. One can employ the in- 
dex theorem for open spaces to estimate the number of these zero energy 
modes [T71IT8] . The index of an operator of the form 



D= fill 




is deflned by index(D) = fc+ — fc_, where fc+ (A;_) is the number of zero modes 
of L [L"^). Its absolute value gives a lower bound to the total number of zero 
modes, fc+ + A;_, oi D. It is possible to bring Hamiltonian fITOl) in the form flTTl) 
by conjugating it with 1 (g) [^(l e^"'^'^^^)], where S{aa (Tb) = CTb® Ca for all 
the Pauli matrices including ctq = 1. Then the index of Hamiltonian flTUj) is 
given by 

index(D) = -^/ rfS-e'^-V^^M!!^, (12) 

which corresponds to the degree of $« if it is considered as a mapping from 
S"^ to S^. For the single Yang-Mills monopole given in (Q with charge 1 we 
have index(D) = 1 indicating that the system has (at least) one zero mode. 



Similarly, it is possible to simulate higher charge monopoles, the explicit forms 
of which are given in [19j . 

It has been demonstrated [20] by analytic calculations that the monopole 
configuration (^ possesses a unique normalized, isolated and non-degenerate 
zero-energy mode. Its wave function is given by '^e=o ~ {4'^'^\ O)'^? where (Z)^"-* 
is a two component spinor 



t\r) 



NieM- f dr'\]-^{r') + A{r')\){5u52n - S2An) (13) 

Jo ^2, J 



with i = 1,2 and A^ a real normalization constant. This solution corre- 
sponds to a wave function that is well localized around the position of the 
monopole, r = 0. The presence of a single zero mode in the spectrum of the 
Dirac operator, together with the symmetry under fermi-number conjugation, 
\^Hc = cTg (g, cr2(c"2)nm^''™''*, of Hamiltoniau fITOl) signals the fractionalization 
of the fermion number [2^ . Thus, the presence of a classical topological defect 
in the configurations of the Yang- Mills and scalar fields causes the fermionic 
field to take fractional quantum numbers. This can be detected in the lat- 
tice fermion system by comparing the fermion densities with and without the 
monopole. A fermion density difference that corresponds to the localized wave 
function (1131) provides a clear signature of the charged fractionalization. 



3 The two dimensional case 



For completeness we would like to present the generation of the Yang-Mills 
field in the two-dimensional case. For that consider a square lattice with tun- 
nelling fermions governed by a similar Hamiltonian as in ([1]). The transition 
couplings are given by X^j = xij + ^Xij > with 

Xi,i+x = ~~7 Xi,i+y — ~'^\^) 

and 



AV AX 

AXM = (-l^f + ^e-(^-^-)-S 
_ A| 

Ay 
Axi,i+x = (-l)-— e-(^+-^-)-', (14) 

where A^ = A^ — iAg , the Fermi points are K+ = (0, 0) and K_ = (0, vr) 
and the Brillouin zone, 5, is given by k^ G [— vr/2, 7r/2] and ky G [— vr, vr]. The 



homogeneous tunnelling couplings Xij give rise to a pair of two-dimensional 
Dirac equation corresponding to the two Fermi points. Under the conjugation 
by f/ = diag(l, a^) the composite Hamiltonian becomes 1 ® 7 ■ p, where 
7 = (0-3, (Ti). The distortion A^i i of the local site energies gives rise to 

From the distortion A^i^i+x+y we obtain 

Ai7(*+^) = -^a3®a3 
while from Axi i+x we have 

It is easy now to show that under conjugation by U the overall Hamiltonian 
produces the two dimensional Dirac operator coupled to an SU{2) Yang-Mills 
field D^ = (1 ® 7 • p - T'^ ® 7 ■ Aa)^. 

It is an exciting possibility that such a lattice Hamiltonian could be gener- 
ated with optical lattices and ultra-cold atom or molecule technology [22j. 
In particular, the generation of the vr-flux through the square plaquettes of 
optical lattices is possible by imprinting phases with laser assisted tunnelling 
processes [23] . The presented method for generating Yang-Mills fields has the 
advantage that it does not employ the internal states of the atoms that de- 
mand additional laser fields [24] and, thus, it requires a significantly simpler 
experimental setting for its physical realization. For example, the distortion 
of the tunnelling couplings is expected to appear naturally in optical lattice 
systems due to amplitude variations of the corresponding laser fields. Alterna- 
tively, controlled spatially profiles of laser fields can be easily produced in the 
laboratory without increasing the complexity of the experiment. Applications 
of non-Abelian fields that facilitate the robust simulation of the quantum Hall 
effect with ultra-cold atoms have been reported in |25j. 



4 Conclusions 



In conclusion, we considered a fermionic system with discretized three spatial 
dimensions described by the Lagrangian 



where time is considered to be a continuous variable. By a particular choice of 
the transition couplings of the fermions we demonstrated that the low energy 
limit of this system is faithfully described by Dirac fermions coupled to a 
triplet scalar field, also known as a Higgs field, and a Yang-Mills SU{2) non- 
Abelian gauge field. These fields can be chosen to encode Yang-Mills-Higgs 
monopole configurations. In particular, a single monopole with charge 1 can 
give rise to an isolated zero-mode. We demonstrated that such a setting gives 
rise to fractionalization of charge in three dimensions. Moreover, we derived 
the analytic relation between the lattice fermions wave function and the Dirac 
ones, thus making it possible to detect the topological effect of the effective 
gauge theory by probing the behavior of the constituent lattice particles. It is 
intriguing to see if such systems in two or three dimensions can be simulated 
in the laboratory. 

The presented simulation opens up numerous possibilities. It would be inter- 
esting to generate kinetic terms for the Yang-Mills field and augment it to 
an independent dynamical entity in a quantum field theory. In this context 
the study of the full symmetry group of the lattice and the emergence of the 
correct massless spectrum in the quantum continuum limit is highly desirable. 
Alternatively, a proof of stability of the quantum order of these non-Abelian 
systems as well as the encoding of chiral fermions would bring this effort closer 
into a full blown simulation of the standard model. 
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